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ABSTRACT 

As the evolutionary link between the radio millisecond pulsars (MSPs) and the low 
mass X-ray binaries or intermediate mass X-ray binaries, the millisecond X-ray pulsars 
(MSXPs) are important objects in testing theories of pulsar formation and evolution. 
In general, neutron stars in MSXPs can form via core collapse supernova (CC channel) 
of massive stars or accretion induced collapse (AIC channel) of an accreting ONeMg 
WD whose mass reaches the Chandrasekhar limit. Here, in addition to CC and 
AIC channels we also consider another channel, i.e., evolution induced collapse (EIC 
channel) of a helium star with mass between 1.4 and 2.5Mq. Using a population 
synthesis code, we have studied MSXPs arising from three different evolutionary 
channels. We find that the Galactic birthrates of transient MSXPs and persistent 
MSXPs are about 0.7—1.4 x 10“'* yr“*. Our population synthesis calculations have 
shown that about 50%—90% of the MSXPs have undergone CC channel, about 10%— 
40% of them have undergone EIC channel, and the MSXPs via AIC channel are the 
least. 
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1 INTRODUCTION 

Millisecond pulsars (MSPs) are a population of old neutron 
stars (NSs) which have been detected in the radio, X-ray and 
gamma ray wavelengths of the electromagnetic spectrum. 
Most of them are in binary systems. Depending on the 
nature of the companion star MSPs are believed to form 
from either low mass X-ray binaries (LMXBs) or interme- 
diate mass X-ray binaries (IMXBs) ( i . e. recycling scenario, 
iBhattacharva fc van den Heuvelill99ll : [Taurisll201ll '). In this 
standard scenario, the NS accrete the matter from the 
companion star, and its spin frequency becomes very high 
on a much longer timescale because of angular momentum 
conservation. During this process we can detect the system 
as an X-ray source and towards the end of this process as 
a millisecond X-ray pulsar (MSXP). MSXPs are divided 
into two subclasses: accretion powered MSXPs and nuclear 
powered MSXPs. The former is an X-ray pulsar powered 
by accreted material releasing gravitational energy, and the 
latter is one powered by accreted material burning in an 
unstable thermonuclear flash. On observations, there are 
type I bursts in nuclear powered MSXPs. However, these two 
classes overlap, burst oscillations were observed also from 
eight accretion powered MSXPs dPapitto et al.|[2014h . Not 
only accretion powered MSXPs but also nuclear powered 
MSXPs must accrete matter from their companions. In our 
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work, considering that both of them include accreting NSs, 
we do not distinguish them, and call accretion and nuclear 
powered MSXPs as MSXPs. Once the companion star can 
not fill its Roche lobe and the mass transfer terminates, the 
NS is observed as a recycled radio MSP. This evolutionary 
scenario has been supported by the very recent discovery of 
”a transition between a rotation powered and an acc retion 
powered state in a binary MSP” JPapitto et ahllioisl l. 

As for 5 Aug 2014, about 200 MSPs in the field and 
120 MSPs in globular cluster were knowr{)|. There ar e 14 
accretion powered MSXPs in iPatruno fc Watts! (l2012l f . 16 
nuclear powered MS XPs in IWatts et al.l ll2008ll . respectively. 
iPapitto et al.l (l2013l i found a MSP (IGR J18245-2452) which 
was transiting from accretion to rotation powered emission. 
This system also includes an accreting MSP. In our work, 
MSXPs are the binary systems including accreting MSPs. 
Therefore, only 31 MSXPs are observed. Figure [T] shows 
the distribution of spin periods (Pg) and orbital periods 
(Porb) of MSXPs and radio binary MSPs. If MSXPs are the 
progenitors of radio binary MSPs, we have to explain the 
followings: (i)Where do the radio binary MSPs with orbital 
periods longer than 100 hours come from? (ii) Where do the 
MSXPs with orbital periods shorter than 1 hour evolve to? 

Usually, the NS is formed via core collapse (CC) 


^ https://apatruno.wordpress.com/about/millisecond-pulsar- 
catalogue/ 



















2 Zhu et al. 





Figure 1. The distributions of spin periods (Pg) and orbital 
periods (Porb) of MSXPs (circles) and radio binary MSPs 
(triples). The observat ional data of th e form er come from 
iPatruno &: Watt3 ||2012| ) and IWatts et al.l ||200S| ) , and the data 
of the latter come from ATNF pulsar catalogue. 


supernova (SN) of a massive star, and it receives additional 
high kick velocity (~ 100 km s“^). Many literatures found 
that the evolutionary scenario from LMXBs or IMXBs to 
MSPs via CC channel cannot explain the origin of binary 
MSPs whose orbital periods are longer t han 10 0 day s 
(e. g., IPfahl et "^120031 : iLiu fc ChenI 12013 ). iMichell (1 19871 ) 
considered that an ONeMg white dwarf (WD) accretes 
sufficient material from its companion, and can form a 
rapidly spinning NS via accretion induced collapse (AIC) 
when its mass reaches the Chandrasekhar limit. Considering 
angular momentum conservation, it is possible that accret¬ 
ing ONe Mg WD directly collapses into M SP via the AIC 
scenario (iBhattacharva fc van den Heuveil Il99ll ). Assumed 
that the AIC scenario r esults in a low kick velocity (« 10 
km s“^), iHurlev et al.l (I2OI0I ) demonstrated that the AIC 
channel cannot be ignored in radio binary MSPs formations 
and that some radio binary MSPs in wide orbits (up to 
abo ut 1000 days) are bes t explained by AIC scenario (Also 
see ISmedlev et al.l l2014l ). IChen et al.f ll201ll) investigated 
whether eccentric binary MSPs can form via AIC channel. 
However, both of them did not consider MSXPs. 

In the above literatures, authors d id not consider tha t 
NS can be formed via the third scenario. iMivaii et al.l lll980l ) 
considered that ONe core of single star can collapse to a 
NS by electron capture on ^'^Ne and ^"^Mg. This process is 
called as ECSN. Because ECSN is largely very dependent 
on stellar composition , mass-loss rate or internal mixing, 
lEldridee fc Tou3 (l2004l) pointed out that the mass range of 
stars under g oing i t is largely uncertain. For solar metallicity, 
IPoelarendsI (l2007l ) suggested that the initial-mass range of 
stars between 9.0 and 9.2Mq can undergo ECSNe . It is 
between ~ 9.8 and 9.9 Modot in iDohertv et akl (l2015l) . For 
a single star, the ECSN channel is limited to a very narrow 
initial mass range. Its contribution to the total gravitational 
coll apse SN is negligible. _ 

IPodsiadlowski et al.l ll2004h discussed the binary evo¬ 
lution examples. The star in a binary system may ex¬ 
perience a process of mass transfer, and it losses whole 
convective envelope prior to central helium ignition and 


becomes a naked helium star. Due to the lack of hydrogen- 
rich convective envelope, the star does not undergo the 
second dredge-up, w hich results in a m ore mass i ve he ¬ 
lium core. Based on Nomotd (Il984l) and iNomotd (Il987l) . 
IPodsiadlowski et al.l (2004) considered that the naked he¬ 
lium st ars with masses between 1 .4 and 2.5 Mq undergo 
ECSN. IPodsiadlowski et al.l ll2004h pointed out that these 
naked helium stars correspond to initial masses of single 
stars about 8 — 11 Mq. Therefore, comparing to a single 
star, the mass range of a star in a binary system for E CSNe 
occurring is large. Following llvanova et al.l (l2008l ). this 
channel is called as the EIC scenario. In their work, NSs in 
binary systems can be formed via three scenarios: CC, AIC 
and EIC. C onsidering dynamica l events (collisions, mergers 
and so on), llvanova et ^ 1I2OO8I) investigated the formation 
and evolution of LMXBs and MSPs in globular clusters, and 
found that AIC and EIC are important formation channels 
for NSs in globular clust ers. Including single and binary 
systems, iKiel et ^ 1I2OO8I) simulated the evolution of the 
Galactic radio pulsars. Compared with the observed diagram 
of spin period vs. spin period derivative for MSPs, they 
found that low kick velocity for AIC and EIC can form a 
greater number of MSPs. Both of them focused on radio 
MSPs population in binary systems. Up to now, there is no 
detailed theoretical work to simulate how MSXPs form via 
the above three scenarios. 

Recently, MSXPs have attracted the attention of as¬ 
tronomers. Many X-ray and gamma ray space missions like 
RXTE, XMM-Newton, INTEGRAL, Chandra, Swift and 
HETE have discovered many observational properties of 
these objects, s uch as the orbital periods, spin period, X- 
ray luminosities (I Watts et al.|[200^: lAlt amirano et al.l I2OI0I : 


|Pa.truno fc WatO 20121 : Papitto et akl 2013h . To our knowl- 

edge, there is no investigation of population synthesis for 
MSXPs to date. As the evolutionary link between radio 
MSPs and NSs in LMXBs/IMXBs, MSXPs have played an 
important role in testing theories of pulsar formation and 
evolution. In addition, the studies of MSXPs will deepen our 
knowledge about the accretion physics and the properties 
of NSs. In this sense, it is necessary to perform a detailed 
population synthesis of MSXPs. 

In this work, we concentrate on MSXPs, and attempt 
to investigate the properties of MSXPs formed via the three 
evolutionary channels. Section 2 presents our assumptions 
and describes some details of the modelling algorithm, and 
Section 3 gives the main results and discussed in detail. Our 
main conclusions are summarized in Section 4. 


2 MODELS 

iHurlev et al.1 1 I 2 OO 2 I) gave the rapid binary star evolution 
(BSE) code which is often u sed in many population synthesis 
codes. [Kiel fc HurlevI 1 I 2 OO 6 I) updated BSE code. In our work, 
we use BSE code in which there are many input parameters. 
If we do not specially m ention any input par ameter, its 
default value is taken from iKiel fc HurlevI (I 2 OO 6 I ). 

The common envelope (CE) evolution and the NS 
formation and evolution are very important factors in our 
simulation. They have great effects on the population of 
MSXPs. The following three subsections give descriptions. 
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2.1 CE Evolution 


In binary systems, CE evolution is very important. Un¬ 
fortunately, the understanding about it is still very poor. 
Usually, people considered that dynamically unstable mass 
transfer results in CE evolution. In BSE code, there are 
two input parameters to model CE evolution. The first 
parameter is the criterion for dynamically unstable mass 
transfer, q^. If the mass ration of the components [q = 
Afdonor/fUgainer) is larger than gc, the binary system ex¬ 
periences CE e volution. Howev e r, the value about q^ is 
still uncertain jHan et al. I l200ll . I2OO2II. Con s iderin g non¬ 
conservative Roche lobe overflow, iHan et (l200 lll foun d 
that qc may be up to about 1.3 (also s eelChen fc Hanll2008ll . 
Very recently, IPavlovskii fc Ivanoval ll2015l l showed that qc 
can varies from 1.5 to 2.2 if the recombination energy in 
the super adiabatic layer of g iant d onor s was conside r ed. In 
addition, iKing fc Beaelm^ (Il999l l and IXauris et al.l (l2000l l 
suggested that the CE evolution is probably avoided in any 
binary in which the convective envelope of the donor star is 
not too deep if accretion rate is super-Eddington. In order to 
discuss the effects of qc on MSXPs po pulation, we s i mulat e 
three cases as the follows: (if Following iHurlev et all ll2002ll . 
we take gc=3 (0.695) if the donor is a main sequence (MS) 
star with mass larger (lower) than 0.7 Mp, respe ctively. If 
the donor is a giant we take qc as ll Webbinkll 1988ll 


qc = 0.362 -I- 


1 

3X (l-Mc/Mdonor)’ 


( 1 ) 


where Me is the core mass of the dono r; _ 

(ii) For giant donor, qc = 1.5 llHan et al.l i200ll '): [Chen fc 
I 2008ll : ?), while other qcS are the same with (i); 

(iii) CE evol ution only occurs when t h e donor in a binar y 
is a giant llKing fc Begelm^ 1 19991 : iTauris et al.l I2OOOII . 
Therefore, for giant donor, gc=Eq(IT]), while in order to avoid 
CE evolution other qcS are equals to 100. 

The second parameter is the efficiency of the binary 
orbital energy used to expel the envelope of donor, ace. 
However, based on a balance eq uation for the orbita l 
angular momentum of binary svstem. fNelemans et al.l (l200Cll l 
introduced another parameter 7 to describe CE evolution. 
The former is called as a-algorithm, and the later is called as 
7-algorithm. Similarly, the value of 7 is also very uncertain. 
In order to discuss the effect of CE evolution on MSXPs 
population, we carry out different simulations in this work. 
In a-algorithm, aceAce = 1.0, and 7 = 1.5 in 7-algorithm. 
Here Ace is a stellar structure parameter, and it depends 
on the evolutionary stage of the donor. As iNelemans et ^ 
(l200(tl suggested, the 7-algorithm can be applied only in 
the first CE evolution for binary systems in which two 
low-mass stars have similar masses. Considering that CE 
evolution is crucial in close binaries, we use 7-algorithm to 
simulate the first CE evolutions for all binary systems. When 
binary systems undergo the second CE evolutions, we use a- 
algorithm. 


2.2 Formation Channels of NSs and Kick Velocity 

As described in the introduction, there a re three channels for 
the formation of NSs in binary systems ll Ivanova et al. I l2008l: 
iKiel et ^l2008l l: (i) CC channel for a star with MS mass 
M/Mq ^11; (ii) AIC channel for an accreting ONeMg WD 


when its mass reaches the Chandrasekhar limit; (iii) EIC 
channel for naked helium stars with masses between 1.4 and 
2.5M©. In BSE code, these naked helium stars originate 
from initially more massive stars in binary systems with 
masses in the range of 8 — 11 M© whi ch is consistent with 
the ra nge of stellar masses quoted by iPodsiadlowski et al.l 
(l2004h . It is question able if an MSP wou ld form directly from 
any kind of channel dTauris et al.ll20l^ . In our simulations 
all NSs formed via three channels need an accretion phase 
and are subsequently spun up to become MSXPs. The 
reasons are discussed in section 3.3.2 

The evolution of accreting ONeMg WD in AIC channel 
is poorly known. For simplicity, we use the same way of 
calculating the evolution of a CO WD by accretion of 
hydrogen-rich material from its companion. The amount 
of matter retained by accreting WD is critical for success 
of AIC channel. For the rete ntion efficiencies for hydrogen 
accumulation (?7h), we follow lHachisu_ef©aT ([199^) (also see 
iHan fc Podsiadlowskil (l2004l l and Lii_et_^ (l2009|)). Based 
on the model of optically thick wind. lKato fc Hachisul ll2004l ~) 
calculated the retention efficiencies for helium accumulation 
(^He). Considering theoretical and observational evidence 
again st existence of optica l ly th ick winds from accreting 
W D dLepo fc van Kerkwiil312013[ ). we used 77He calculated 
by IPiersanti et al. ( 20141 ') who took into account existence 


of Roche lobes in binaries instead of assuming optically 
thick winds. However, in order to discuss the effects of the 
retention efficiencies on MSXPs population, we also used 
Kato-Hachisu ret ention efficiencie s in a case. In addition, 
from the work of Ijose fc Hernand (Il998l l. it is known that 
ONe WDs have about 2 times higher retention efficiencies. 
Therefore, we carried out two simulations in which the 
retention efficiencies of accr eting ONe WD s were 2 or 1/2 
times of those calculated bv IPiersanti et al] ||2014| ~). 

The limits of stable accumulating matter also determine 
the mass growth of accreting WDs. The limit s for accreting 
CO WDs depend on the WDs masses (e. g. , iHachisu et al.l 
I1999I : iKato fc Hachisul I2OO4I : IPiersanti et al.ll2014l 'l. and the 
larger is CO WDs mass, the higher is the limit. Usually, 
compared to CO WDs, ONe WDs are massive. Therefore, 
in our work, we extrapolated the limits for CO WDs to ONe 
WDs although they may be difference. 

Not only observations but also theoretical simulations, 
it is very difficult to giv e a certain value for kick velocity 
received by nascent NS. iPfahl et al.l 1I2OO2II suggested that 
the nascent NSs via CC channel have high kick velocity (~ 
100km s“^), while the NSs born in EIC and AIC have low 
kick velocity(-^ 10 km s“^). Like in our study abou t the 
donors of NSs in persistent LMXBs dZhu et al.ll2012l l. the 
distribution of kick velocity (iik) is given by 



iHobbs et al.l (l2005l l investigated the proper motion of 233 
pulsars, and found that a Maxwellian distribution with (Jk = 
265 km s“^ can describe well the velocity distribution of 
these pulsars. However, in order to check the effect of Uk on 
MSXPs pop ulation, we take Uk = 190 and 400 km s“^ for 
CC channel llHansen fc Phinn^ll997l L while = 20 and 
50 km s”*^ for EIC and AIC channels in different simulations, 
respectively. 
















































































































4 Zhu et al. 


2.3 Evolutions of NSs Spin and Magnetic Field 


The evolutions of NSs spin and magnetic field are mainly 
determined by the interaction of the NS with surrounding 
materials. This interaction has become the most important 
but still is an open problem. For the ev olution of NSs 
spin, w e use a convenient way elaborated bv lLimmov et all 
(Il992tl . Spin evolution (spin-up or spin-down) of a NS in a 
binary system can be conveniently described by an angular 
momentum conservation equation 


dloj 

dt 


- .^sd 5 


( 3 ) 


where 7 is NS momentum of inertia, Tfsu and K^d are spin- 
up and spin-down torques, respectively. We treat the NS’s 
spin evolutions in the same way as in our p revious study 
for symbiotic X-ray binaries dLii et al.l l2012l l. The values 
of spin-up and spin-down torques exerted on a rotating 
magne tized NS are summarized in Table 2 of iLii et al] 
(l2012tl . In this work, the value of the index n, which reflects 
different rotational distribut ions of matter inside t he shell, is 
taken 2 (|Shakura et al.l[2012h ffor further details see lLii et al.l 

iI2oI3) ). 

There is no an exact mech anism to describ e the decay of 
NSs magnetic field. Following iLii et aP l|2012h . we assume, 
for an accreting NS, the decay of magnetic field depends 
exponentially on the am ount of accreted mat ter. We use 
the formulae suggested bv lOslowski et all (1201 ll i. For a non¬ 
accreting NS, we assum e that the decay of magnetic fie ld re¬ 
sults from Ohmic decay ([Goldreich fc Reise neggeill l992l'). an d 
the decay formulae is given by Eg. (2) in 1^1 et_alj _l[200^ 
Follow ing iFaucher-Giguere fc Kasol (12006 ) and lPopov et al] 
(l2010li . a log-normal distribution of the initial magnetic 
fields of nascent NSs is assumed. The median value of the 
distribution is 12 and its dispersion equals 1. 

In LMXBs/IMXBs, it is possible that not all matter 
from their companions via Roche lobe flows or stellar winds 
is accreted by NSs. We introduce a parameter /3 which is 
the fraction of transferred matter accreted by the NS, and 
the rest of the transferred matter is lost from binary system. 
The lost matter takes away the specific angular momentum 
of t he gainer. The value of / 3 has been usually set to 
0.5 I Podsiadlowski et al.lll9^ . l2002l : iNelson fc RappaportI 
I2OO3 ). In our work, in order to compare the influence of /?, 
we set P = 0.25 and 1.0 in different simulations. In many 
systems, the mass transfer rates may exceed the Eddington 
limit. However, not all of materials are accreted. In this 
work, we assume that a star can accrete material up to the 
Eddington limit, and any excess material is lost from the 
binary system. 


2.4 Persistent and Transient X-ray Sonrces 

On observation, the MSXPs can be classified into 
persistent and transient X-ray sources. In all 31 
known MSXPs, t here are 24 tran si ent s o urces and 7 
persi st ent sourcesdPatruno fc Wattj I2OI2I : IWatts et al.l 
l2008l : IPanitto et ahl I2OI3I I. lii this paper, we consider 
the thermal disk instability and make a distinction 
between persistent and transient sources. The thermal disk 
instability depends on the mass-accretion rate of NS. The 
X-ray binary is a transient source when the mass-accretion 
rate is lower than a certain critical value. Merit. That is. 


if Mns < Merit in Roche overflow accretion, the system is 
a transient source, here Mns is the mass-accretion rate of 
NS. If Mns > Merit or wind-fed accretion, the sy s tem i s 
a persistent source. According to iBelczvnski et al.1 (l2008ll . 
for the hyd rogen-rich an d heav ier-el ement disks we us e 
the works of IPubus et al.1 (llQflflI l and iMenou et al.1 (l2002l l , 
respectively: 


' 1.5 X 10i®M-°-‘77^'iC'j-° ® g/s, H rich 
5.9 X 10i®M-° «^77^ ®2ag:f g/s, He rich 
1.2 X 10i®M-° '^^7?l-2Vg:f g/s, C rich 
5.0 X 10i®M-° ®®7?^ o®a§:f g/s, O rich 


where Mns is the NS mass in solar mass, Rd is the maximum 
disk radius (it is approximately equal to 2/3 of accretor 
Roche lobe radius) in 10^® cm, constant Ci = C/(5 x 10““*) 
in which C is a radiation parameter and its typical value 
is 5 X 10““^, ao.i = a/O.l in which a = 0.1 is a viscosity 
parameter. 


3 POPULATION OF MSXPS 

Using the method of population synthesis for binary popula¬ 
tion, several initial input parameters (initial mass function, 
initial mass ration of binary components and initial binary 
separation) must be given. For the primaries in binary 
sy stems, their ma s ses ar e given by the initial mass function 
of iMiller fc Scal3 (Il979tl . For initial mass ration of binary 
components, we adopt a flat d istribution jKraicheva^£_^ 

I I989I: [Goldberg fc Mazehlll994l l. Following lYungelson et ah 
(1993), we assume that initial binary separation, Ao, is 
between 10 and 10® T?©, and has a fiat distribution over 
log Ao. 

In this work, we considered eight cases in which dif¬ 
ferent input parameters are listed in Table [T] In order to 
investigate the birthrates of MSXPs, we assume simply a 
constant star formation rate over last 13 Gyr and that 
there is one binary with Mi ^ O.SMq born per year 
in the Galaxy J^en fc TutukovlflOSl : lYungelson et al.lll993l : 
iHan et ^ ]fi9^ In every case, we simulate the evolutions of 
2 X 10® binary systems. This gives a statistical error between 
about 1% and 2% for the birthrate of the Galactic MSXPs 
via different channels. 


3.1 Birthrates and Average Lifetime 

The birthrates and the average lifetime of the transient 
MSXPs and the persistent MSXPs with different kinds of 
donors in the Galaxy are summarized in Tables [2] and [3l 
respectively. 

In our simulations, the total birthrate of the transient 
MSXPs arising from the three channels is from 6.9 xl0~® 
(Gase 2) to 1.4 xl0“"^yr“^ (Case 3). The total birthrate of 
the persistent MSXPs is from 8.4 xl0“® (Case 4) to 1.3 
xl0~"^yr~^ (Case 1). Considering CC and AIC channels, 
iHurlev et al.1 (I2OI0I I studied radio binary MSPs, and showed 
that their birthrate is between 5.5 xl0~‘*yr“^ and 2.2 
xl0“®yr“^ in the simulations with different input param¬ 
eters. According to the theories of radio MSPs formation 
and evolution, the radio emission is activated when the 
companion of NS can not fill up its Roche lobe and the mass 
transfer terminates. That is, the MSXPs will eventually 
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Table 2. Different models of the transient MSXPs’ population. The first column gives the model number according to Tabled 
Columns 2 to 7 give the birthrates and the average lifetime of MSXPs with different kinds of donors, respectively. Total birthrate 
and lifetime are showed in columns 8 and 9, respectively. NS + MS means that accreting NS has a main sequence donor, NS + RG 
means that accreting NS has a red giant donor, NS + WD represents that accreting NS has a WD donor, and ’—’ means that we 
do not get any sample in 2 x 10® binary systems. 



NS+MS 

NS+RG 

NS+WD 

Total 


Cases 

Birthrate 

Lifetime 

Birthrate 

Lifetime 

Birthrate 

Lifetime 

Birthrate 

Lifetime 


(yr”^) 

(yr) 

(yf^) 

(yr) 

(yr”^) 

(yr) 

(yr”^) 

(yr) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 





cc 





Cases 1, 6—8 

3.5x10-® 

2.5x10® 

1 .6x10-® 

7.0x10® 

1.3x10-® 

1 .2x10® 

6.4x10-® 

2.3x10® 

Case 2 

1.7x10-® 

2.7x10® 

2.5x10-® 

7.3x10® 

1 .6x10-® 

1.9x10® 

2 .1x10-® 

2 .6x10® 

Case 3 

5.9x10-® 

3.3x10® 

9.4x10-® 

3.1x10® 

5.2x10-® 

2 .1x10® 

7.4x10-® 

3.2x10® 

Case 4 

2 .0x10-® 

2.9x10® 

1.7x10-® 

1 .0x10® 

9.0x10-® 

1 .0x10® 

4.6x10-® 

2 .8x10® 

Case 5 

4.4x10-® 

2.5x10® 

2 .0x10-® 

5.3x10® 

2 .1x10-® 

9.0x10’ 

8.5x10-® 

2.3x10® 





AIC 





Case 1 

— 

— 

l.OxlO-’’ 

8 .1x10® 

7.2x10-® 

8 .1x10® 

7.4x10-® 

8 .1x10® 

Case 2 

— 

— 

O.OxlO-’’ 

4.5x10® 

1.4x10-® 

3.8x10® 

1.4x10-® 

3.8x10® 

Case 3 

— 

— 

1.3x10-® 

2 .0x10’ 

— 

— 

1.3x10-® 

2 .0x10’ 

Case 4 

— 

— 

7.6x10-® 

1 .6x10® 

4.5x10-® 

1.4x10® 

1 .2x10-® 

1.4x10® 

Case 5 

— 

— 

7.3x10-’’ 

8.5x10® 

1.3x10-® 

8 .0x10’ 

1.4x10-® 

8 .0x10’ 

Case 6 

4.8x10-'^ 

4.3x10® 

5.4x10-’’ 

9.0x10’ 

1.4x10-® 

5.8x10® 

1.5x10-® 

1.3x10® 

Case 7 

1 .6x10-'^ 

3.6x10® 

4.4x10-’’ 

1.7x10’ 

1.4x10-® 

3.8x10® 

1.5x10-® 

6.9x10® 

Case 8 

— 

— 

1 .0x10-’ 

3.2x10® 

1.4x10-® 

3.8x10® 

1.5x10-® 

3.8x10® 





EIC 





Case 1, 6—8 

1 .8x10-® 

1.3x10® 

9.0x10-® 

6 .2x10® 

1 .2x10-® 

4.4x10® 

3.9x10-® 

1 .1x10® 

Case 2 

2.5x10-® 

1.3x10® 

3.7x10-® 

1.4x10® 

4.9x10-® 

4.5x10® 

3.4x10-® 

1.3x10® 

Case 3 

2 .8x10-® 

1 .0x10® 

1 .6x10-® 

2.9x10® 

4.7x10-® 

6 .1x10® 

6 .0x10-® 

6 .0x10® 

Case 4 

1 .8x10-® 

1 .0x10® 

8.7x10-® 

4.4x10® 

1 .2x10-® 

3.6x10® 

3.9x10-® 

8 .1x10® 

Case 5 

1.9x10-® 

1 .0x10® 

1 .2x10-® 

4.6x10® 

1 .6x10-® 

3.9x10® 

4.7x10-® 

7.8x10® 


Table 3. Similar to Table[2] but for the persistent MSXPs’ population. 



NS+MS 

NS+RG 


NS+WD 

Total 


Cases 

Birthrate 

Lifetime 

Birthrate 

Lifetime 

Birthrate 

Lifetime 

Birthrate 

Lifetime 


(yr"^) 

(yr) 

(yr"^) 

(yr) 


(yr-i) 

(yr) 

(yr"^) 

(yr) 

(1) 

(2) 

(3) 

(4) 

(5) 


(6) 

(7) 

(8) 

(9) 






CC 





Case 1, 6^—^8 

2 .8x10-® 

4.5x10® 

4.0x10-® 

1 .1x10® 


1 .6x10-® 

1 .2x10® 

4.8x10-® 

4.3x10® 

Case 2 

1 .6x10-® 

5.1x10® 

1 .0x10-® 

1 .1x10® 


2.7x10-® 

1.3x10® 

2 .0x10-® 

5.0x10® 

Case 3 

4.5x10-® 

4.8x10® 

5.5x10-® 

5.3x10® 


8.7x10-® 

1.3x10® 

5.9x10-® 

4.7x10® 

Case 4 

1.5x10-® 

5.7x10® 

8 .2x10-® 

1.5x10® 


2 .8x10-® 

1.3x10® 

5.1x10-® 

5.0x10® 

Case 5 

3.4x10-® 

4.7x10® 

7.1x10-® 

9.7x10® 


1 .6x10-® 

1 .2x10® 

5.7x10-® 

4.5x10® 






AIC 





Case 1 

— 

— 

1.5x10-’ 

4.8x10® 


4.5x10-® 

1 .2x10® 

4.6x10-® 

1 .6x10® 

Case 2 

— 

— 

3.1x10-’ 

2.7x10® 


8.5x10-® 

1 .2x10® 

8 .8x10-® 

1.3x10® 

Case 3 

— 

— 

1.3x10-® 

3.1x10® 


— 

— 

1.3x10-® 

3.1x10® 

Case 4 

— 

— 

1.7x10-® 

1.4x10® 


2.3x10-® 

1 .2x10® 

4.6x10-® 

1.3x10® 

Case 5 

— 

— 

7.0x10-’ 

4.7x10® 


6.7x10-® 

1 .2x10® 

7.4x10-® 

2 .2x10® 

Case 6 

4.8x10-’ 

3.7x10® 

4.9x10-’ 

6 .6x10® 


8 .6x10-® 

1 .2x10® 

9.6x10-® 

2 .2x10® 

Case 7 

1 .6x10-’ 

3.7x10® 

4.2x10-’ 

5.9x10® 


8.5x10-® 

1 .2x10® 

9.1x10-® 

1.3x10® 

Case 8 

— 

— 

1 .0x10-’ 

3.9x10® 


8.5x10-’ 

1 .2x10® 

9.5x10-’ 

1 .6x10® 






EIC 





Case 1, 6^—^8 

5.6x10-® 

2.3x10® 

1 .2x10-® 

2 .8x10® 


1 .1x10-® 

1.9x10® 

7.9x10-® 

2 .2x10® 

Case 2 

8.9x10-® 

1 .6x10® 

3.6x10-® 

2 .2x10® 


4.1x10-® 

2 .1x10® 

9.8x10-® 

1 .6x10® 

Case 3 

5.6x10-® 

4.1x10® 

1 .8x10-® 

2 .8x10® 


3.8x10-® 

1.5x10® 

6 .2x10-® 

2.9x10® 

Case 4 

5.6x10-® 

2.3x10® 

1 .2x10-® 

2 .8x10® 


1 .1x10-® 

1.9x10® 

2.9x10-® 

1.7x10® 

Case 5 

5.8x10-® 

2.3x10® 

2 .0x10-® 

2 .1x10® 


1.4x10-® 

1.9x10® 

4.0x10-® 

1 .6x10® 
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Table 1. Parameters of the models for MSXPs’ populations. The 
first column gives model number. Column 2 shows the algorithm 
of CE evolution. Column 3 gives the dispersion of kick-velocity 
distribution for CC channel, and column 4 for AIC and EIC 
channels. The critical mass ratio for dynamically unstable mass 
transfer, is given in column 5, in which (i), (ii) and (iii) 
are described in subsection 2.1. Column 6 shows the retention 
efficiencies for heliu m accumulation (rjHp) , in which P and K- 


^He are 

calculated bvlPiersanti 
respectively. 

et al.l 1 I 2014 I') 

andiKato 

& Hachisul 

Cases 

CE 

C7k(km s 

o-J(kmi 

3“^) <?C 

7He 

Case 1 

^ace'^ce ~ 1-0 

190 

20 

(i) 

P-7He 

Case 2 

^ace-^ce — I'O 

400 

50 

(i) 

P-7He 

Case 3 

7 = 1.5 

190 

20 

(i) 

P-7He 

Case 4 

^ce-^ce — 1-0 

190 

20 

(ii) 

P-7He 

Case 5 

^ace-^ce — I'O 

190 

20 

(iii) 

P-7He 

Case 6 

^ace-^ce — I'O 

190 

20 

(i) 

K-7He 

Case 7 

<aceAce = 1-0 

190 

20 

(i) 

2xP-7He 

Case 8 

<aceAce — I'O 

190 

20 

(i) 

ixP-7He 


evolve to the radio MSPs, which means that the birthrate 
of the radio MSPs is about —10 ~^ yr~^ in the Galaxy . 

Our results consist with those in iHurlev et al.l ll2010f) . 
However, the birthrate of the radio MSPs establis hed from 
their radio properties, is about 3—5xl0~^ vr~^ ([Lorimej 
I 2 OO 5 I : iFerrario fc Wickramasingh3 l2007l : IStorv et all 2007ll , 
a factor 10—100 times lower than our calculations. This 
means that we encounter a problem of overproduction of 
MSXPs. 

The overproduction problem may results from the 
followings: 

(i) Selection effects. The sample of observed MSXPs is 
heavily biased towards the brighter objects that are the 
easiest to detect. The observed MSXPs are only the tip of 
the iceberg of a much larger underlying population. 

(ii) The calculation of mass-transfer rates. The lifetimes of 
the MSXPs depend on mass-transfer rates. In our work, 
we use BSE code to calculate the mass-transfer rates. 
Compared to full evolutionary computations, BSE code may 
give inaccu r ate r esults. Using BSE code and MESA code, 
IChen et al.l (l2014l l calculated the mass-transfer rates of some 
binary systems, respectively. They found that BSE code 
choices lower mass-transfer rates than that in MESA code 
when the donors are red giants. Our work may overestimate 
the lifetimes of MSXPs. 

(iii) The X-ray active lifetimes of the transient MSXPs. 

In general, the transient MSXPs are relatively faint X- 
ray sources. Their typical quiescent luminosities are about 
10 ^ 2 — j^q 33 g-i peak X-ra y luminosities during 

outbursts are ab out 10^® erg s“^ (IPapitto et al.l I 2 OI 1 I : 
iPaizis et al. I [112). During the quiescent phase, it is very 
difficult to observe the transient MSXPs as X-ray sources. 
The duty cycles (the fraction of time during the outburst) 
are very uncer tain. Empirically the d uty cycles is not 
exceed ing 1 % jBelczvnski et all l2008l i. Ivan Haaften et ^ 
(l 2012 tl used theoretical approach to calculate duty cycles 
which depend on the mass, radius and mass-accretion 
rate of the gainer and binary separation. They found 
that the duty cycles in some ultra-compact X-r ay binaries 
decrease below 0.1%. However. lYan &: Yul ll2015li performed 
a statistical study of the outburst properties of 110 bright 


X-ray outbursts in 36 LMXB transients. They found that 
the duty cycles of these transients with NSs are about 1%— 
10%, and the average duty cycles are 3.5%. 

A similar problem appears in the study of LMXBs: 
about 1000 —10000 strong X-ray systems ar e predicted 
by population synthesis codesJPfahl et aL||2003^. while less 
than 200 LMXBs are observed ( Liu et al. 2007ll . There¬ 
fore, in order to estimate their number in the Galaxy, 
we assume that the observed pr obability of a t ransient 
MSXPs approximately equal 3.5% llYan fc Yull2015h . Then, 
we estimate that there are about 3700 (case 2) — 9800 (case 
3) observable transient MSXPs in the Galaxy. According to 
our calculations, about 53% (Gase 2)— 86 % (Gase 3) of the 
transient MSXPs arise from GC channel, and about 14% 
(Case 3)—42% (Case 2) of them arise from EIC channel, 
while the transient MSXPs arising from AIC channel are less 
than 3% and can be negligible because of too short lifetime. 

Similarly, based on Table |31 there are about 340 (case 
2) — 500 (case 3) persistent MSXPs in the Galaxy. About 
30% (Case 2)—70% (Case 4) of the persistent MSXPs arise 
from CC channel, about 10% (Case 4)—40% (Case 2) of 
them arise from EIC channel and about 10% (Case 3)—30% 
(Case 2) of them arise from AIC channel. In our work, there 
should be 4000—10000 MSXPs in the Galaxy. However, up 
to now, there are only 31 known MSXPs in the Galaxy. 
Like previous studies on LMXBs, we also overestimate the 
number of MSXPs. 

In our work, an important result must be noted. 

Altho ugh it is unclear whether all MSXPs become 

MSPs dTaurisI l 2012 ll . we still assume that all MSXPs 
should eventually evolve to radio MSPs. Taking Case 1 
as an example, we find that the birthrate of the radio 
MSPs formed v ia CC channel is 6. 4 x 10“® yr“^ which is 
consistent with iHurlev et al.l (I 2 OIC 1 I I . The birthrates of the 
radio MSPs via AIC and EIC channels are 7.4 x 10“® and 
3.9 X 10“® respectively. Usually, all nascent radio 

MSPs have low magnetic fiel ds. Their lifetimes as radio 
sourc es approximately equal llEerrario fc Wickramasinghd 
l2007li . This means that about 58% of radio MSPs result 
from CC channel, 7% of radio MSPs result from AIC 
channel, and 35% of radio MSPs result from EIC channel. 
Therefore the EIC channel is one of the most important 
channels in the formation of the radio MSPs, and cannot 
be ignored in studies of the MSPs. 


3.2 Effects of Parameters 

In our work, different input parameters have effects on 
the populations of MSXPs. Compared with Case 1, the 
parameters ak and cr^ are increased from 190 to 400 km s“^ 
and from 20 to 50 km s“^ in Case 2, respectively. The larger 
ak and al are, the more difficultly a binary survives after 
SN. Therefore, the birthrate of MSXPs in Case 2 decreases. 

In general, OceAce = 1.0 in a-algorithm means that 
a binary orbital period after CE evolution should shorten 
to 1% —10% of the orbital period at beginning of Roche 
lobe overflow. While, 7 = 1.5 in 7 -algorithm means that 
a binary orbital period is approximately constant during 
CE evolution. In our simulations, different algorithms of CE 
evolutions have different effects on the MSXPs via different 
channels. For CC and EIC channels, the birthrate of the 
MSXPs in Case 3 is larger than that of in Case 1, however. 
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for AIC channel, the result is opposite. The reason is that, 
for AIC channel, an ONeMg WD should accrete sufficient 
matter from its companion star before it collapses to a 
rapidly spinning NS. The distances of post CE systems 
with the a-algorithm are narrower than those with the 7 - 
algorithm, which is favor for ONeMg WDs accreting more 
matter. So, for AIC channel, the birthrate of the MSXPs 
in Case 1 (a-algorithm) is larger than that of in Case 3 ( 7 - 
algorithm). 

Compared to Case 1, Qc is larger for giant donors in 
Case 4. It means that CE occurs more difficultly in the 
latter. Therefore, there are more NS + RG systems and 
less NS + WD systems in Case 4. The assumption of 
in Case 5 mainly avoids CE evolution when donor is on 
Hertzsprung Gap, and the material of the donor’s envelope 
transfers to gainer. When the donor evolves into giant, its 
envelope becomes CE. Based CE evolutionary theory, less is 
the mass of CE, more can binary survive after CE formation. 
In Case 5, there are more MSXPs. According to Cases 1, 
4 and 5, Qc affects the birthrates and lifetimes of MSXPs 
within a factor of 1.5. 

The MSXPs’ birthrate via AIC channel mainly depends 
on the parameter ^He (retention efficiency for helium accu¬ 
mulation ) in this work. Compa ring Cases 1, 6 and 7, we 
find that the rjne c alculated bv iK ato fc Hachisul (120041 ') is 
higher than that bv IPiersanti et al.l ( 2014|l . and the former 
is approximately two times as high as the latter. From 
? 7 He = 2 X P-rjiie in Case 7 to r^He = 1/2 x P-rjue in Case 8 , 
the birthrate of MSXPs via AIC channel reduces by about 
10 times. 


3.3 Properties of MSXPs via the Three 
Evolutionary Scenarios 

The present work focuses on the MSXPs populations via 
CC, AIC and EIC scenarios. In this subsection we discuss the 
properties of MSXPs via the three evolutionary scenarios. In 
the following we take Case 1 as an example to discuss some 
properties (orbital periods, spin periods, donor’s masses, 
mass-accretion rates) of MSXPs populations. 


3.3.1 Initial Orbital Period and Donor’s Mass: vs. 

Mi 

Figure [ 2 ] shows the distributions of the initial orbital periods 
Porb vs. the initial masses of secondaries (They will become 
the donors of MSXPs) for the progenitors of MSXPs. 
Compared with the progenitors of transient sources, the pro¬ 
genitors of persistent sources have larger initial secondaries’ 
masses but they have similar orbital periods. On average, 
among the three evolutionary scenarios, the progenitors of 
MSXPs via EIC scenario have the shortest initial orbital 
periods because they must undergo CE evolution to form 
naked helium stars when the primaries just leave MS stage, 
while the progenitors of MSXPs via AIC scenario have the 
longest initial orbital periods because the primaries must 
have enough space to form ONe WDs. 

Most of binary systems in Figure [2] undergo violent 
binary interaction, and then evolve into MSXPs. Compared 
with their progenitors, MSXPs have very shorter orbital 
periods which are showed in Figure O The left graph is 


for the transient MSXPs, and the right is for the persistent 
MSXPs. For the transient sources, there are three different 
zones for Porbs’ distributions. The left zone (Porb < 1 hour) 
presents the ultra-compact MSXPs, that is, the donors in 
these MSXPs are WDs or naked helium stars. The transient 
MSXPs with MS donors mainly are in the middle zone 
(1 < Porb < 50 hours), and the right zone (Porb > 50 
hour) is mainly composed of the transient MSXPs with 
red giant donors. All of three types of donors appear in 
the transient MSXPs arising three channels. Obviously, the 
orbital periods of transient MSXPs via AIC channel are 
shortest because in our simulatious they undergo twice 
stable mass transfers: one is that the accreting ONe WD 
collapses into NS, another is that the accreting NS is 
spin-up into MSP. Therefore, in our work, AIC channel 
is no contribution to MSXPs with long orbital periods. 
iHurlev et al.l (l201Cllf assumed that accreting ONe WD can 
directly collapse into MSP without a spin-up process. In 
their work, the contribution of binary MSPs via AIC channel 
is significant. 


3.3.2 Orbital Period: Porb 

Similarly, as the right graph in Figure [3] and Table |3] show, 
the most of the persistent MSXPs have MS donors, and the 
ultra-compact persistent MSXPs only is about 2% of the 
whole persistent MSXPs. Comparing the left graph with the 
right one in Figure O the orbital periods of the transient 
MSXPs can be wider than those of the persistent MSXPs. 
The main reason comes from our definition for the transient 
and the persistent MSXPs (See Eq. ([4|)). 

In all known MSXPs, the histogram of their orbital 
periods observed are plotted in the bottom panels of Figure 
[31 As the middle and the bottom panels in Figure [3] show, 
our results cover all observational samples well. However, 
two aspects in our results are worth noting: 

(i) In all known MSXPs, the orbital period (18.95 hours) 
of Aql X-1 is the most long. However, according to Figure 
[T] many radio binary MSPs have orbital periods longer 
than 100 days. In fact, in our simulation, there are some 
MSXPs with orbital periods longer than 100 days. However, 
the duration of they staying in X-ray stage is too short. 
Therefore, they are negligible and hardly are showed in 
the histogram. Figure [4] shows the distributions of orbital 
periods of MSXPs when they terminate mass transfer 
and become radio binary MSPs. If all radio binary MSPs 
originate from MSXPs, it is difficult to use our model to 
explain the formation of r adio binary MSPs w ith long orbital 
period. Similarly done in iHurlev et al.l (l201(]f) . we also do a 
simulation in which we assume that radio binary MSPs can 
be directly formed via AIC channel without spin-up process. 
Our result is plotted by dashed line in the middle panel of 
Figure [4] Obviously, our results still do not explain why so 
many radio binary MSPs have orbital periods longer than 
1000 hours. 

(ii) Based on Tables [2] [3] and Figure [S] the ultra-compact 
transient MSXPs only is about 4% of the whole transient 
MSXPs, and the ultra-compact persistent MSXPs only 
is about 2% of the persistent MSXPs. However, in the 
known 24 transient MSXPs, there are 5 ultra-compact 
transient MSXPs, and in the known 7 persistent MSXPs, 
there is 1 ultra-compact persistent MSXP. Obviously, we 
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Figure 2. Gray-scale maps of the distributions of the initial orbital periods P^rh vs. the initial masses of accreting NSs’ companions in 
the MSXPs. The left is for the transient sources, and the right is for the persistent sources. 


underestimate the ratio of the ultra-compact MSXPs to the 
total MSXPs. The possible causes will be discussed in the 
next subsection. 

In additi on, where do thes e ultra -compact MSXPs 
evolve into? IPelove fc Bildst^ ll2003ll considered the 
effect of ideal gas pressure, degeneracy pressure and 
Coulomb attraction on the WD radius, and found that 
some ultra-compact MSXPs can evolve into radio binary 
MSPs with orbital periods 70 — 90 minutes within 5 
— 10 Gyr. Recently, having thought the propeller effect 

and disk instabi lity _ in ultra-compact X-ray binaries, 

Ivan Haaften et ^ (|2012h considered that they evolved 
into low mass ratio binaries with orbital periods 70 — 
80 minutes. In the ultra-compact MSXPs, we use a zero- 
temperature WD mass-radius relation, and orbital changes 
mainly originate from gravitational radiation and mass 
variations. In our work, with the donor’s mass decreases 
which results in the decrease of mass transfer rate, the 
ultra-compact MSXPs evolve into the faint ultra-compact 
X-ray binaries, and hardly evolve into radio binary MSPs 
in Hubble time. Therefore, as Figure [4] shows, there is no 
radio binary MSP with orbital period shorter than 1 hour 
in our work. However, if we assume that MSPs can directly 
form via AIC channel without spin-up process, there are 
some radio binary MSPs with orbital periods shorter than 
1 hour (See the dashed line in the middle panel of Figure 
SJ. They may evolve into the ultra-compact MSXPs due to 
gravitational radiation within 10 Gyr. Therefore, we should 
observe some radio binary MSPs with orbital periods 
shorter than 1 hour. However, to our knowledge, there is 
no known radio binary MSP with orbital period shorter 
than 1 hour. We considered that there may be two kinds of 
explanation: 

(i) If ONe WD can directly collapse into a MSP via 
AIC channel, the newborn MSP may have high magnetic 
field because no a c creted material buries magnetic field. 
Ivan Paradiis et ^ l|l997t l considered that the bursting 
pulsar GRO J1744-28 originates from a massive ONe 
WD via AIC (Also see IXu fc Lil l2009h . Although GRO 
J1744-28 is not a MSP (its pulsar period is 467 ms), 
Ivan Paradiis etliL 1 119971) assumed that it was a MSP when 


it was born. According to the observati o nal d ate of pulse 
period and spin-up rate, iFinger et al.l lll996ll estimated 
that the magnetic field of GRO J1744-28 is of order 10^^ 
G. Therefore, the newborn MSP via AIC channel may 
be a magnetic field of ~ 10^^ G. Treating the pulsar as 
a rotating magnetic dipole, the surface magnetic field 
strength B — 3.2 x 10^®(PP)^^^ G, and the characteristic 
age Tc = P/(2P). Then, we can estimate that Tc ~ 10® yr 
if the newborn MSP via AIC channel has a spin period of 
about 10 ms. Compared to normal MSPs with low magnetic 
field (P ~ 10®“® G) and long characteristic age(rc ~ 10® 
yr), the lifetime of the newborn MSP via AIC channel as 
a radio MSP is very short. Therefore, it is difficult to find 
them. 

(ii) The process of WD collapsing into a MSP via AIC 
channel may be very complex, and the angular momentum 
of ONe WD can be lost by some mechanism. For example, 
the SN via AIC scenario ejects some matter which can carry 
out more angular momentum. However, to our knowledge, 
there is not any simulation or observation to refer it. 

Therefore, it is possible that NS via AIC channel can 
not directly become MSPs, or it is difficult to observed these 
MSPs even though they can be directly formed via AIC 
channel. However, if the NSs via AIC channel undergo spin- 
up processes to become MSPs, they should be similar with 
other MSPs. Considering the above reasons, we assume that 
the MSPs via AIC channel must undergo spin-up processes 
in the present work. 


3.3.3 Orbital Period and Donor’s Mass: Porb vs. M 2 

Figure [5] gives the distributions of the donors’ masses vs. 
orbital periods of the MSXPs. Compared with Figure (2] the 
progenitors of MSXPs have undergone violent mass transfer. 
The distribution of the companion stars’ masses for AIC 
channel is narrower than those of for CC and EIC channels. 
The interpretation for this distinction is that AIC binary 
systems may have experienced significant accretion phases, 
and a significant amount of material has been transferred 
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log ^orb(l^OUrs) 

Figure 3. The distributions of the orbital perio ds of MSXPs . The l eft is for the trans i ent sources, and the right is for the persistent 
sources. The data of observed MSXPs come from lWatts et aTI l|2008lb [Patruno &: Watt3 l|2Q12li : [Papitto et al.l l|^13l 'l. 
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Figure 4. The distributions of the orbital periods of MSXPs 
when they terminate mass transfer and become radio binary 
MSPs in Hubble time. The dashed line in the middle panel 
represents radio binary MSPs directly formed via AIC channel 
without spin-up process. The orbital periods observed in all 
known radio binary MSPs are plotted in the bottom panel. The 
data of observed MSXPs come from ATNF pulsar catalogue. 

to NSs, then the masses of companion stars are usually 
relatively small. In ultra-compact MSXPs, the distribution 
of donors’ masses are cut into two regions. In the upper 
region (M 2 > OAMq), the most of donors are naked helium 
stars , while they are WDs in the down region. For the 
MSXPs with orbital periods longer than 100 hours, donors 
mainly are red giant and they have low masses (« 0.5Mq) 
in order to avoid dynamical mass transfer. 

3 . 3.4 Spin Period: Ps 

The spin periods are one of the most important physical 
parameters for MSXPs. Figure [6] shows the distributions of 
the spin periods of NSs in MSXPs. For the transient MSXPs 
(See the left graph in Figure [6]), there are two peaks. The 
left peak around about 1.6 ms originates from the transient 


MSXPs with MS or red giant donors, and the right peak 
around about 20 ms from the transient MSXPs with WD 
donors. Usually, the mass transfer between two degenerate 
binaries mainly is driven by gravitational radiation, and the 
mass transfer between a degenerate star and a MS or red 
giant mainly is driven by stellar evolution dHurlev et al.l 
l2002li . In our work, the time of the former occurring is 
much longer than that of the later occurring, which means 
that the NSs in transient MSXPs with WD donors have 
lower magnetic fields because of the decay of magnetic fields 
dLu et al.ll2012l '). A lower magnetic field results in a shorter 
spin period. For the distribution of spin periods of NSs in 
the persistent MSXPs (See the right graph in Figure [B)), 
there is a peak at Pg ~ 5.6 ms. As Figure [6] shows, the 
majority of NSs in the transient MSXPs have spin periods 
between about 1.4 and 6.0 ms, while the spin periods of NSs 
in the persistent MSXPs mainly are between about 6.0 and 
20 ms. In our simulations, the transient MSXPs usually have 
shorter spin periods than those of the persistent MSXPs. 
The main reason is that, compared with the NSs in the 
persistent MSXPs, the NSs in the transient MSXPs have 
longer accretion histories which results in lower magnetic 
fields. As the bottom panels in Figure [6] show, our results 
are in agreement with the observational samples. 

3.4 Corbet Diagram of MSXPs 

According to the distributi on of spin periods Ps vs. orbital 
periods P...h. fCorbed (Il986l l divided high-mass X-ray bina¬ 
ries into three classes. Therefore, Corbet diagram is very 
important for us to know X-ray binaries. In the known 
31 MSXPs, there are 18 transient sources and 4 persistent 
sources whose orbital periods and spin periods are measured. 
They are plotted in Figure [7| in which the distributions 
of the spin periods Ps vs. the orbital periods Porb for the 
MSXPs simulated in our work are given in grey-scale. As 
Figure 0 shows, whether observed or simulated MSXPs are 
cut into two zones: ultra-compact and normal sources. For 
the transient and persistent MSXPs, the observed normal 
sources are covered by the simulated results in all cases. We 

















































































10 Zhu et al. 



-1.0 0.0 1.0 2.0 3,0 4.0-1.0 0.0 1.0 2.0 3.0 4.0 -1.0 0.0 1,0 2.0 3.0 4.0-1.0 0.0 1.0 2.0 3.0 4.0 
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Figure 5. Similar with Figure [2] but for grey-scale maps of the distributions of the orbital periods Paih vs. the donor’s masses in the 
MSXPs. The left is for the transient sources, and the right is for the persistent sources. 




-3 - 2.5 -2 - 1.5 -3 - 2.5 -2 - 1.5 

log P,is) log P„(s) 

Figure 6. The distributions of the spin periods of NSs in the MSXPs . The le ft is for the trans i ent sources, and the right is for the 
persistent sources. The data of observed MSXPs come from lWatts et ahl ll2n08l 'l : fPatruno fc Watt3 ll2012h : rPapitto et al.l ll201.'il) . 


can not suggest which one is better. However, all simulations 
underestimate the rotation speeds of the NSs in ultra¬ 
compact MSXPs. As the last section mentioned, in this 
work, the mass transfer between two degenerate binaries 
mainly is driven by gravitational radiation. However, it 
is very difficult to accurately simulate the mass transfer 
between two degenerate binaries. Our models may underesti¬ 
mate the mass-transfer rate between WDs and NSs. A high 
mass-transfer rate is favor for forming MSXPs. Therefore, 
underestimating the mass-transfer rate between WDs and 
NSs is just the reason that we underestimate the ratio of 
the ultra-compact MSXPs to total MSXPs. 


4 CONCLUSIONS 

MSXPs are thought to be the progenitors of radio MSPs. 
As the evolutionary link between radio MSPs and NSs in 
LMXBs/IMXBs MSXPs have played an important role in 
testing theories of pulsar formation and evolution. In this 


paper, we have investigated in detail the MSXPs arising 
from three different evolutionary channels (including the 
EIC channel, which has almost been ignored in previous 
studies of radio MSPs), using the population synthesis 
method to the binary evolutions. 

Our simulations predict that the Galactic birthrates of 
the transient MSXPs and the persistent MSXPs are about 
0.7—1.4 X 10“"^yr“^. The observable number of the transient 
and the persistent MSXPs in the Galaxy are about 4000— 
10000 and 300—500, respectively. We have discussed the 
properties of the MSXPs, such as orbital period, spin period, 
masses of the companion stars, and mass-accretion rates 
of the NSs. Our results about the distributions of spin 
periods and orbital periods for MSXPs are similar to the 
observations. 

Our population synthesis calculations have shown that 
about 53%—86% of the transient MSXPs have undergone 
GC channel, about 14%—42% of them have undergone EIC 
channel, and the transient MSXPs via AIG channel are 
negligible. The proportions of the persistent MSXPs via the 
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Figure 7. Similar with Figure!^ but for grey-scale maps of the distributions of the spin periods Ps vs. the orbital p eriods Pprh for 
MSXPs. The left is for the transient sources , and the right is for the persistent sources. The data of observed MSXPs ll Watts et al.|[2003: 
IPatruno Watt3l2012l: IPapitto et al■|[2013^ are plotted by pentagram. 


above three channels are abont 30%—70%, 10%—40% and 
10%—30%, respectively. According to the theories of MSPs 
formation and evolution, the MSXPs will eventually evolve 
to the radio MSPs. Our simulations show that about 58% 
of radio MSPs result from CC channel, 7% of radio MSPs 
result from AIC channel, and 35% of radio MSPs result from 
EIC channel. Therefore the ETC channel is one of the most 
important channels in the formation of MSXPs and radio 
MSPs, and cannot be ignored in studies of MSPs. 


However, the investigation of MSXPs population in¬ 
volves many uncertainties (See section 2). Although we 
try our best to choose commonly accepted paradigms, our 
results can not explain observations in many respects, such 
as: an order of magnitude mismatch between observed and 
predicted theoretically numbers of MXBPs, and underesti¬ 
mating the spin periods of accreting NSs in ultra-compact 
MSXPs and their numbers. The possible reasons are as 
the follows: (i) erroneous treatment of binary evolution, 
including CE evolution, accreting NS evolution, and so 
on; (ii) incorrect lifetimes involving mass-transfer model 
and uncertain duty cycles of transient MSXPs; (iii)not 
considering the observational selection effects. In a short, 
there still is a long way to go for completely understanding 
MSXPs. 
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